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Abstract 

Gravitational fluctuations along the line-of-sight from the surface 
of last scatter to the observer distort the microwave background in 
several related ways: The fluctuations deflect the photon path (grav- 
itational lensing), the decay of the gravitational potential produces 
additional fluctuations (ISW effect) and scattering off of hot gas in 
clusters produce additional fluctuations (Sunyaev-Zel'dovich effect). 
Even if the initial fluctuations generated at the surface of last scatter 
were Gaussian, the combination of these effects produce non-Gaussian 
features in the microwave sky. We discuss the microwave bispectrum 
as a tool for measuring a studying this signal. For MAP, we estimate 
that these measurements will enable us to determine the fraction of 
ionized gas and to probe the time evolution of the gravitational po- 
tential. 



1 Introduction 

Over the past few years, cosmologists have emphasized that measurements of 
the two-point correlation function of CMB anisotropics can yield a wealth of 
information about physical conditions in the early universe. In this article, 
we will discuss how non-linear physics can produce a detectable non-Gaussian 
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signature in the microwave background. We concentrate on cases in which 
gravitational lensing between here and the surface of last scatter couple via a 
Limber's equation to produce a non-Gaussian signal. The amplitude of this 
signal is sensitive to the evolution of gravitational potential fluctuations at 
low redshift. Thus, measurements of the non-Gaussianity will complement 
measurements of the two point correlation function by providing information 
about physical conditions in the low-redshift universe. 

In §2, we discuss cross-correlating the effects of gravitational lensing and 
the Integrated Sachs- Wolfe (ISW) effect and cross-correlating the effects of 
gravitational lensing with the Sunyaev-Zel'dovich (SZ) effects. Because there 
are equal numbers of hot and cold spots, gravitational lensing alone does 
not produce a detectable three-point signal [[rj. However, the combination 
of lensing with other low redshift effects will produce an observable non- 
Gaussian signature. In §3, we compute the signal/noise for this effect and 
discuss these cosmological applications. A companion paper provides a 
more general discussion of the bispectrum, including derivation of the angle- 
averaged bispectrum, a method for measuring the bispectrum and estimating 
the signal to noise, and a fully worked example of the bispectrum produced 
by the Rees-Sciama effect. 



2 Computing the Bispectrum 

The angle-averaged bispectrum is the spherical harmonic transform of the 
temperature three point function, 



Bhhh = E Li „ 3 / didriuin(r(i)T(m)T(n)> 

x ^: mi (i)^ m2 (m)^ ms (n) , (1) 



where T(n) is defined as the temperature anisotropy (with mean zero) at a 
particular point in the sky. We expand the temperature out into three terms: 

T(n) = aimYi m (n) + ]T a im V0(n) ■ VY lm {h) + T, r (n) (2) 

Im Im 

where the first term is the CMB fluctuations from the surface of last scatter, 
the second term describes the effects of lensing and the third term describes 
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the effects of secondary processes at low redshift. We shall assume throughout 
this paper that the third term is first order (e.g. oc <&o(k)), but clearly, second 
order secondary processes would also give rise to an observable bispectrum, as 
is the case with the Rees-Sciama effect considered in the companion paper [0]. 
For this article, we will consider two sources of secondary anisotropy (both 
of which are first order), the ISW effect and the Sunyaev- Zel'dovich effect. 
Combining equations (|l|) and (||) yields, 

x £ y Im ft(a lm o? m ,>(Ve(m).V17 m ,(m)T fc .(n)> (3) 

plus an additional 5 terms representing the various permutations of the ob- 
servation angles. 

We can immediately integrate over 1: 



Bhhh = E t 2 ) jd^YU^U^) (4) 



h h h 
rrii ni2 rris 

x Ql (Ve(m).V^ mi (m)T ir (n)) 



We may note that by integrating by parts: 

|rfmV 2 ^ m3 (m)6(m)^ mi (m) = - J rfmV^ m2 (m) • Ve(m)^ mi (m) 

- J dmVr i2 * m2 (m) • Vl£ mi (m)e(mX5) 

where the surface integrals drop out, since we are integrating over a complete 
unit sphere. Combining this expression, with similar expressions involving 
the Laplacian of 6(m) and YJ* (m), we find: 

Bhhh = ~ T I h h h ) f dhY*(n) Cll (6) 

x 



J dm([v 2 r i2 * m2 (m)0(m)^ mi (m) 



- F4Jm)V 2 0(m)^ mi (m)-r Z2 * m2 (m)e(m)V^: mi (m^ T lr (n) 
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Since G(m) and T} r (n) are scalar quantities, their expectation value will 
depend only on their spatial separation: 

(e(m)TUn)) = E(^)Wm-fi) (7) 

lm 

Inserting equation (Q) into equation (0) and integrating yields, 



k(k + 1) - h{h + 1) - Z 3 (fe + 1) (2/i + 1)(2Z 2 + 1)(2Z 3 + 1) 
Bhhh ~ 2 V 4^ 

x(o o) C ^ (8) 
Plus 5 additional terms, reflecting the various permutations of l\, I2 and I3. 

2.1 Cross-correlating Lensing with the ISW Effect 

If the universe is not flat with Q m = 1, then at late times, gravitational po- 
tential fluctuations decay This decay generates additional microwave back- 
ground fluctuations: 

T ISW {h) = 2 / dr(p(T,nr) 
Jo 

= 2/ dr— -3$ o (k)exp(zk-nr)0(r) (9) 

JO (,27TJ 

where the gravitational potential fluctuations have been expanded out in 
Fourier space, n is the direction of photon propagation, and r is conformal 
(comoving) lookback time. If a line-of-sight passes primarily through high 
density regions, then it will be on average hotter than lines-of-sight that pass 
primarily through low density regions. 

Foreground fluctuations in the gravitational potential also distort the 
photon geodesies, 

0(m) = 2 f° dr^^ V±0(r, mr) (10) 
Jo r 

dr 3-/c$o(k)(k — (k ■ m)rh) exp(ik • mr)0(r) 

(27r) to 

= v±e(m) 



4 



This effect is surprisingly large; a typical photon's trajectory has been de- 
flected roughly several arcminutes from its initial path. 

By measuring the three-point function, we can cross-correlate these two 
effects, 

(T ISW (h)Q(ra)) = - f k 2 dkP(k) H dr H dr>(r)^^ 

7T J JO JO TTq 

x0(r') Y,Ji(kT) 3l (kT')Y* n (*h)Y lm (h) , (11) 

l,m 

where we have used the Rayleigh expansion of the Fourier mode into spherical 
harmonics: 

e 1 ^ = A^ l 3i(kr)Y; m (v)Y lm ii) ■ (12) 

lm 

Since most of the contribution to the integral comes from when r ~ r' , 
we approximate <j) by its value at that time and then integrate over r'\ 

p (L _|_ I) 

b\ sw = -_-A^Z f P{k)kdk f dr<P{r) 3l {kr)^r) T -^- (13) 

In our analysis, we compute this integral numerically up to I — 50. However, 
for large /, we can evaluate this integral quickly by noting that ji(kr) ~ 
5{kr — l)Jir/2l, yielding: 

bj SW ~ -8/ dT<J)(T)<P(T)^-^P(l/T) (14) 

To illustrate the behavior of this coefficient, Figure [I] shows a plot of b\ versus 
I for an Q m = 0.3, A = 0.7 cosmology. 

2.2 Cross-correlating Lensing with the Sunyaev-Zel'dovich 
Effect 

Filaments, groups, and clusters will distort the microwave background through 
both the SZ effect and gravitational lensingQ. The amplitude of the effect 
depends upon the integrated pressure fluctuations along the line-of- sight, 

T sz (n) = [ adr5p e (r) (15) 

m P c z J 
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where ot is the Thompson cross-section, m e is the electron mass and Sp e is 
the fluctuations in electron pressure. 

The electron pressure is the product of the electron density and the elec- 
tron temperature. On large scales, we parameterize this as Sp e = b gas pA, 
where b gas is the gas bias factor (which is expected to be ~ 4 as dense regions 
are hotter than low density regions Q), p is the mean gas pressure and A is 
the fluctuation in the dark matter density. This approximation only holds 
on scales in which linear theory holds. However, as we will show, the SZ- 
lensing bispectrum probes scales of k~ x > 100 Mpc. Numerical simulations 
(R. Cen, private communication) find that the density weighted temperature, 
kT e ~ 1/(1 + z) keV for z < 3. Since the gas density falls as (1 + z)~ 3 , the 
pressure drops roughly as (1 + z)~ 2 . Based on this, we approximate the large 
scale SZ fluctuations as, 

T*'W = -2^ (16) 
J T a 

where yo = aTn e0 kT e oT /m e c 2 . For a flat universe with H = 65km/s, and 
^ionized = 0.05 T e0 = 1 keV, y Q = 1.2 x 10 -5 . This approximation is likely to 
be valid on scales larger than the non-linear scale. 

We can relate the dark matter density to the potential and rewrite equa- 
tion (|T§, 

and now compute the SZ — lensing cross-correlation term: 

(T 5Z (n)6*(m)> = 

x 

As I gets large, the spherical Bessel functions can be treated as approximating 
delta functions such that: 

Ukr') ^ 5(1 - kr')^^^ (19) 



dkk 4 P(k) / dTdT(j)(T)(f)(T' 



Y,Ukr)ji(kr')Y; n Xni)Y lm (n) 

lm 



T 



TTn 



(18) 
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Making this approximation for r' and integrating from — oo to oo yields: 



This may be further simplified by making the same approximation for 
ji(kr) and integrating over k: 



Note that for very small values of r (low redshift), P{l/r) oc r~ 7 , and the ker- 
nel of the above equation goes to 0. It is only in the large wavelength regime 
in which P(k) oc k~ 5 where the integral is maximized. This occurs around 
A = 100/i _1 Mpc for most cosmological models, and thus, the assumption 
that pressure and density perturbations are proportional and evolve linearly 
are well motivated. 

As Figure [I] shows, the Sunyaev-Zel'dovich effect dominates over the ISW 
effect at large I. Since the SZ effect is enhanced by non-linear evolution, 
the estimated signal in equation (pip is likely an underestimate, numerical 
simulations, complemented by improved analytical techniques, will be needed 
to more accurately compute the signal. 

It may be thought that bi will merely change normalization with a change 
of Q m , and thus, yob gas (which clearly only affects the normalization) and fl m 
would be indistinguishable. Since there is little variation in q for variations of 
fl m , it might be thought that the bispectrum shape will be roughly constant, 
and thus, only the normalization could be extracted. In Figure 0, we show 
that this is not the case. For variations in Q m , we see that bi varies not only 
in normalization, but in shape as well, and thus, both yob gas and Q m will be 
able to be extracted almost independently. 

3 Discussion 

Will we be able to detect either effect? We can estimate this by determining 
whether in a Q m = 0.3, Q\ = 0.7 universe, we can use the MAP data to 
reject the hypothesis that there is no non-Gaussianity in the CMB maps: 




(21) 



x 2 = E 



h,h,h 




(22) 
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where we restrict the sum to l\ < I2 < h- For our standard set of parameters, 
X 2 — 9 with most of the signal coming from the SZ effect. Thus, measuring 
the bispectrum will enable us to measure O, ionized and probe the time evolution 
of the gravitational potential. 

In Figure |3|, we plot the value of x 2 for the ISW-lensing bispectrum for flat 
models with varying Q m and As shown in the Figure, the PLANCK || 
satellite will do an excellent job differentiating between the null hypothesis 
and a given value of fl m , giving a 4 a signal for Q m = 1.0. The MAP satellite, 
however, will be unable to differentiate between the two. 

Likewise, this method will be able to potentially probe different equations 
of state. While the Cosmological constant behaves as a component with 
w = P/p = — 1 (our fiducial model), different equations of state produce 
dramatically different rates of potential growth, and hence, we may measure 
the x 2 difference between test values of w and w = — 1 , as shown in Figure [|. 
Note that both MAP and PLANCK provide constraints on w at the 1 a level; 
PLANCK gives w < -0.8, and MAP gives w < -0.4. 

As seen in Figure |3], MAP will be sensitive to this effect as well, and the 
4 year results will yield 1 a uncertainties of 5Q m = ±0.4. PLANCK can 
distinguish SZ from CMB fluctuations, so it will be very sensitive to this 
effect as well. However, calculating its signal to noise is beyond the scope of 
this paper. 

Figure [| shows the additional x 2 f° r a mode with a given value of Z3 (by 
convention, the largest index). Note that the maximal contribution occurs 
near l 3 = 900 for PLANCK. Though the SZ-lensing effect is generally at 
higher mode numbers (see the comparison between the signals as will be 
detectable by MAP), it will nevertheless produce a much stronger signal in 
MAP than will the ISW-lensing effect, since at high I (but less than 600), 
the SZ-lensing coefficient (bi) dwarfs the coefficient for the ISW-lensing effect 
(see Figure p. 

In this article, we estimated the non-Gaussian signature produced by the 
combined effects of gravitational lensing and the ISW effect and the combined 
effects of gravitational lensing and the SZ effect. For most cosmological 
models, the SZ signature is stronger and will swamp the ISW effect signal. 
Detecting the SZ-lensing cross-correlation will determine the mean density 
of ionized gas in the today's universe. Thus, MAP should be able to detect 
the "missing baryons" || |7| . 

By cross-correlating the MAP data with observations of large scale struc- 
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ture and X-ray maps, we can look for additional signatures of low redshift 
physics. Previous papers have discussed cross-correlating the galaxy distribu- 
tion in the Sloan Survey with the temperature maps to look for gravitational 
lensing effect s|| and for the Sunyaev-Zel'dovich produced by both clusters 
and superclustersQ. Boughn et al. [pT0|| have looked for cross-correlations be- 



tween the cosmic X-ray and microwave backgrounds in an attempt to detect 
the late-time ISW effect. Combining these different cross-correlations will 
enable cosmologists to determine the relationship between the mass distri- 
bution and the distribution of gas and galaxies. These observations will also 
probe the time evolution of the gravitational potential and provide a new 
tool for measuring the basic cosmological parameters. 
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Figure 1: The bi coefficients as a function of I. The solid line shows the 
coefficient from the IS W effect , while the the dashed line shows the coefficient 
arising from the Sunyaev Zel'dovich effect. Each of the plots is for an Q m = 
0.3, = 0.7, h = 0.65 cosmology. We can estimate the contribution that the 
different coefficients will have through dimensional analysis. Since q oc l~ 2 , 
the Wigner 3-j symbols oc l~ 1 / 2 (equation H) an d the signal per l 3 is expected 
to go as -B/jZaig x l~ 4 (see equation |22p the x 2 /AZ 3 will be roughly constant if 
bi oc /~ 3 . Thus we have normalized the kernels with an I 3 prefactor in order 
to estimate their importance. 
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Figure 2: As in the previous figure, but the coefficients, bi, are shown only for 
the SZ-lensing coupling. Here we have varied fl m to illustrate that the shape 
of the coefficients, as well as the normalization, vary with Q m , and thus, Q m 
cannot be considered a degenerate parameter with the product, yob gas . 
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Figure 3: The x 2 differences in the ISW-lensing bispectrum, the SZ-lensing 
bispectrum, and their sums between a fiducial model Q m = 0.3, f2 A = 0.7 
and a test model with a flat cosmology and varying Q m . The solid line shows 
the x 2 f° r the PLANCK experiment, while the dashed line shows x 2 for MAP 
4 year results and the the dotted line shows the MAP 1 year result, with 1 
a uncertainty. 
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w =( P a/Pa) 



Figure 4: The same as the previous figure, but with variations in the equa- 
tion, w. For a cosmological constant, w — — 1 (our fiducial model), while 
w = —1/3 is a curvature-like term. Both MAP and PLANCK will be able 
to distinguish between extreme cases; MAP would predict w < —0.4 and 
PLANCK would give w < -0.8. 
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Figure 5: The \ 2 P er mode (l 3 ) as given by the bispectrum. The solid 
line shows the results from the ISW-lensing effect measured by PLANCK, 
while the dashed line shows the SZ-lensing effect as measured by MAP (4 yr. 
results). In general, the signal is at higher wavenumber (smaller scale) for 
the SZ-lensing effect than for the ISW-lensing effect. However, PLANCK's 
increased sensitivity at higher wavenumber distorts this relation in this plot. 
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